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Abstract: Knowledge of soil carbon (C) distribution and its relationship with the environment can 
improve our understanding of its biogeochemical cycling and help to establish sound regional models of 
C cycling. However, such knowledge is limited in environments with complex landscape configurations. In 
this study, we investigated the vertical distribution and storage of soil organic carbon (SOC) and soil 
inorganic carbon (SIC) in the 10 representative landscapes (alpine meadow, subalpine shrub and meadow, 
mountain grassland, mountain forest, typical steppe, desert steppe, Hexi Corridor oases cropland, Ruoshui 
River delta desert, Alxa Gobi desert, and sandy desert) with contrasting bioclimatic regimes in the Heihe 
River Basin, Northwest China. We also measured the 9'Sr/936Sr ratio in soil carbonate to understand the 
sources of SIC because the ratio can be used as a proxy in calculating the contribution of pedogenic 
inorganic carbon (PIC) to total SIC. Our results showed that SOC contents generally decreased with 
increasing soil depth in all landscapes, while SIC contents exhibited more complicated variations along soil 
profiles in relation to pedogenic processes and parent materials at the various landscapes. There were 
significant differences of C stocks in the top meter among different landscapes, with SOC storage ranging 
from 0.82 kg C/m? in sandy desert to 50.48 Ко C/m? in mountain forest and SIC storage ranging from 
0.19 ke C/m? in alpine meadow to 21.91 kg C/m? in desert steppe. SIC contributed more than 75% of 
total C pool when SOC storage was lower than 10 kg C/m?, and the proportion of PIC to SIC was greater 
than 70% as calculated from Sr isotopic ratio, suggesting the critical role of PIC in the C budget of this 
region. The considerable variations of SOC and SIC in different landscapes were attributed to different 
pedogenic environments resulted from contrasting climatic regimes, parent materials and vegetation types. 
This study provides an evidence for a general trade-off pattern between SOC and SIC, showing the 
compensatory effects of environmental conditions (especially climate) on SOC and SIC formation in 
these landscapes. This is largely attributed to the fact that the overall decrease in temperature and increase 
in precipitation from arid deserts to alpine mountains simultaneously facilitate the accumulation of SOC 
and depletion of SIC. 
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1 Introduction 


Soil represents the largest carbon (C) pool (2500 Pg) in the terrestrial ecosystem, storing 
approximately twice the quantity of C in the atmosphere (760 Pg) and vegetation (560 Pg) 
(Schlesinger, 1997; Eswaran et al., 2000; Lal, 2004a). Thus, an increase in C storage capacity of 
soils would be an ideal option for reducing atmospheric CO» (Lal, 2004b; Trumbore and 
Czimczik, 2008). Soil C pool consists of soil organic carbon (SOC) and soil inorganic carbon 
(SIC). The abundance of SOC is often controlled by the balance of natural and anthropogenic C 
inputs (e.g., plant residues and manuring) and outputs through soil organic matter (SOM) 
decomposition (Wynn et al., 2006). The estimated global SOC storage in the top meter ranges 
from 1200 to 1600 Pg (Eswaran et al., 1995; Batjes, 1996). SIC occurs primarily as carbonate 
minerals (Batjes, 1996), such as calcite (CaCO3) and dolomite (CaMg(CO3)2). Theoretically, the 
SIC pool can be divided into lithogenic inorganic C (LIC) and pedogenic inorganic C (PIC). The 
former is inherited from the soil parent material with no chemical change, and the latter is formed 
during the precipitation of Ca?* and/or Mg?* and bicarbonate (West et al., 1988; Lal, 2004c; Liu et 
al., 2014). Consequently, PIC formation might sequester atmospheric CO?» in soils, depending on 
the sources of Са or Mg?* (Monger and Martinez-Rios, 2000). For example, if the Ca?* ог Mg?* 
originates from the dissolution of carbonate minerals during pedogenesis, thus there is no net 
change in C storage (Monger et al., 2015; Schlesinger, 2017); however, if the Ca?* or Mg”* results 
from the weathering of Ca or Mg-bearing silicate minerals, then there will be a net increase in the 
PIC pool (Monger and Martinez-Rios, 2000; Chang et al., 2012; Monger et al., 2015). In addition, 
exogenous sources of Ca?* and/or Mg?* imported from dust, rain and irrigation can also result in a 
SIC sink (Chang et al., 2012; Bughio et al., 2016). The estimated global SIC storage in the top 
meter varies from 700 to 1700 Pg (Eswaran et al., 1995, 2000; Batjes, 1996). The large variations 
in the estimation of global SOC and SIC stocks by different researchers are partly attributed to the 
uncertainties in calculating C storage in regions with contrasting landscapes and bioclimatic 
conditions. Thus, a better knowledge of SOC and SIC distribution and their influencing factors in 
regions with different landscapes and bioclimatic conditions are not only meaningful in 
understanding their regional heterogeneity, but also important for accurately estimating the 
terrestrial C pool. 

Covering 47.2% of the terrestrial area, drylands play an important role in the global C cycle 
(Lal, 2004c; Li et al., 2015). Dryland soils contain 241 Pg of SOC (Eswaran et al., 2000), and 
store at least as much SIC as SOC pool (Eswaran et al., 2000; Lal, 2004c). Previous studies 
suggested that dryland ecosystems have a strong C sequestration potential (Lal, 2004c; Li and 
Shao, 2014; Li et al., 2015). Yet little information is available on the distribution and storage of 
SOC and SIC in the dryland regions, due to their highly heterogeneous environmental conditions. 
The Heihe River Basin (HRB) is the second largest inland river basin in the arid and semi-arid 
regions of Northwest China (Kang et al., 2007; Cheng et al., 2014). Alternating mountains, oases, 
and deserts are typical landscapes in the HRB, showing evident vertical zonal characteristics of 
climate, vegetation and soils (Kang et al., 2007; Yang et al., 2017). So HRB is very suitable for 
studying the C distribution patterns across different landscapes from the alpine to the Gobi 
ecosystems. Yang et al. (2017) found that SOM accumulation, carbonate leaching and 
calcification are the predominant pedogenic responses to soil forming processes in the HRB. As a 
result, the gained, lost and existing SOC and SIC pools as well as the distributions in soil profiles 
across different landscapes have important ecological implications in the inland river basins. 

There have been many studies in the past several decades investigating the storage and 
distribution of SOC at different spatial scales (Batjes, 1996; Schlesinger, 1997; Eswaran et al., 
2000; Jobbágy and Jackson, 2000; Wu et al., 2003a; Lal, 2004a; Wang et al., 2004; Baumann et 
al., 2009; Yang et al., 2016) and their evolution under the influences of natural and anthropogenic 
perturbations in Northwest China (Wang et al., 2004; Wang et al., 2005; Kang et al., 2007; Li and 
Shao, 2014; Lü et al., 2014; Su et al., 2015; Chen et al., 2016; Yang et al., 2016). However, little 
attention has been paid to SIC, in spite of its relative abundance in the arid and semi-arid regions. 
Several studies have demonstrated that the SIC pool in the northern China is one to nine times 
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higher than the SOC stock (Mi et al., 2008; Wu et al., 2009; Wang et al., 2010; Wang et al., 2015; 
Zhang et al., 2015). Often, environments with adequate precipitation favoring SOC accumulation 
have little SIC, while SIC-rich soils occur in the arid and semi-arid regions with lower biomass 
production and little SOC. This suggests that there may be an offset effect between SOC and SIC. 
Thus, studies focusing only on SOC dynamics without consideration of SIC may underestimate 
the C sequestration potential of the terrestrial ecosystem, especially in the arid and semi-arid 
regions where SIC stock is often much higher than SOC stock (Eswaran et al., 2000; Zhang et al., 
2015). 

The objectives of the present study were to investigate the vertical distribution and storage of 
SOC and SIC across different landscapes in the HRB and to identify the relative roles of climatic, 
edaphic and biotic factors in determining SOC and SIC storage and their vertical distribution 
patterns. Results of this study are essential to improve our understanding of regional C cycle in 
relation to the high spatial heterogeneity of climatic conditions, vegetation types and soil 
properties in the arid inland river basins of Northwest China. 


2 Materials and methods 


2.1 Study area 


The Heihe River Basin (HRB) is located in the arid and semi-arid regions of Northwest China 
(37°72'—42°68'N, 97904'—102900"Е), covering an area of approximately 1.28x10? km? (Cheng et 
al., 2014). Geographically, the Неће River can be divided into the upper, middle and lower 
reaches by the Yingluoxia and Zhengyixia hydrological stations, respectively (Fig. 1). The Qilian 
Mountains attach to the upstream areas, the middle Hexi Corridor belongs to the midstream 
regions and the northern Alxa High-Plain belongs to the downstream areas. The Qilian Mountains 
at the upstream areas have elevations ranging from 2000 to 5500 m. The mean annual 
precipitation (MAP) increases from 100 to 500 mm, while the annual mean temperature (AMT) 
decreases from 6?C to —8?C with increasing elevation. The upstream areas can be divided into 
glacier, snow zone, permafrost zone, mountain forest and grassland zone, and steppe zone 
according to the elevation and local topography (Kang et al., 2007; Qin et al., 2013; Yin et al., 
2015). The middle Hexi Corridor with alluvial fans and floodplains is located between the Qilian 
Mountains and the Beishan Mountains. The precipitation is scarce in this region, and the limited 
water resources are collectively utilized in the artificial oases by the highly developed irrigation 
system (Wang et al., 2007; Cheng et al., 2014), where irrigated agriculture is very well developed. 
The lower reaches of the Heihe River flow across the Alxa High-Plain with a mean elevation of 
1000 m. This region is extremely arid with MAP less than 50 mm and annual potential 
evaporation greater than 3500 mm (Yin et al., 2015), where desert and Gobi prevail. 


2.2 Field investigations 


A total of 10 representative landscapes were selected in the HRB (Table 1): alpine meadow (AM), 
subalpine shrub and meadow (SSM), mountain grassland (MG), mountain forest (MF), typical 
steppe (TS), desert steppe (DS), Hexi Corridor oases cropland (HCO), Ruoshui River delta desert 
(RDD), Alxa Gobi desert (AGD), and sandy desert (SD). Field soil survey and sampling were 
conducted during the summer months (July to August) of 2012 and 2013. Three sampling sites 
were set in each landscape, with a total of 30 sites in the whole study area (Fig. 1). Geographic 
coordinates and elevations of each sampling site were recorded using a GPS unit (Juno SB, 
Trimble, USA), and the plant species were also recorded. In each sampling site, a soil pit was hand 
dug to a depth of 120 cm (or bedrock). Soil samples were collected according to genetic horizons 
(Appendix, Table S1). The samples were air-dried and then sieved through a 2-mm nylon mesh. A 
subsample of the air-dried sample was ground using a pestle and mortar to pass through a 0.25-mm 
nylon sieve before laboratory analyses of SOC and CaCO; contents. In addition, undisturbed soil 
samples from each genetic horizon were collected in triplicate using a cutting ring (5 cm in 
diameter and 5 cm in height) for bulk density (BD) determination. The volume percentage of 
coarse fragments (diameter >2 mm) was determined by visual estimation in the field. 
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Fig. 1 Overview of the Heihe River Basin (HRB) and locations of the sampling sites 


Table 1 Characteristics of the 10 representative landscapes in the Heihe River Basin (HRB) 


Altitude МАР AMT А Canopy Parent ; 
Landscape o Dominant species cover А Soil type 
(m) (mm) (°С) (95) material 
3899 324 —4.5 ; z - 75 е 
АМ (4275) (421) (20) Kobresia humilis, Carex tristachya (49) EOD, TIL Histic Cryosols 
3508 332 205 Salix cupularis, Potentilla fruticosa, 57 Cambic 
SMI (+191) (+18) (1.2) Caragana jubata, Kobresia humilis (+13) IOE, SAL Phaeozems 
3212 376 –0.4 Роа annua, Elymus dahuricus, Kobresia 85 Calcic 
MG (+174) (+28)  (=0.9) humilis, Polygonum viviparum (+6) LOE, ALL Kastanozems 
3016 308 0.4 Я ТЕ : T 90 Chernic 
MF (435) (418) (203) Picea crassifolia, Kobresia humili (45) LOE, SAL Phacorens 
2848 272 0.9 Achnatherum splendens, Stipa krylovii, Poa 44 : : 
TS (+118) (+28) (+07) еар (46) LOE Cambic Calcisols 
2039 144 4.9 Sympegma regelii, Pennisetum 13 ; у 
D» (+150) (+12) (40.7) centrasiaticum, Achnatherum splendens (+4) SE Cambio Galutsols 
HCO Les HN 10 Zea mays, Triticum aestivum D HTM Irragric Anthrosols 
(+137) (+21) (+0.6) ' (+5) 
1125 52 8.7 А ; 
RDD (474) (£18) (201) - ALL Skeletic Fluvisols 
997 32 8.7 ; 
AGD (434) G7) (40.1) - EOD, ALL  Fluvic Cambisols 
SD Tey 20 Rc - EOS Calcaric Arenosols 


(+266) (+30) (20.9) 


Note: AM, alpine meadow; SSM, subalpine shrub and meadow; MG, mountain grassland; MF, mountain forest; TS, typical steppe; DS, 


desert steppe; HCO, Hexi Corridor oases cropland; RDD, Ruoshui River delta desert; AGD, Alxa Gobi desert; SD, sandy desert. It 
should be noted that RDD, AGD and SD are barren lands without vegetation. -, no data; MAP, mean annual precipitation; AMT, annual 
mean temperature. Soil types were classified according to the World Reference Base for Soil Resources (IUSS Working Group, 2014). 
Parent material terminology follows the methods of Schoeneberger et al. (2012). EOD, aeolian deposit; TIL, till; LOE, loess; ALL, 
alluvium; SAL, slope alluvium; HTM, human-transported materials (Siltigation in Chinese Soil Taxonomy); EOS, aeolian sands. The 
values in brackets are the standard errors. 


2.3 Soil laboratory analysis 


Soil pH was measured using a digital pH meter in a 1:2.5 soil:water solution. Soil BD was 
measured by drying the undisturbed soil samples at 105?C to a constant weight. Soil particle size 
fractions were analyzed by a laser diffraction particle size analyzer (LS230, Beckman Coulter, 
USA). Prior to the analysis, about 0.2 g soil samples (diameter «2 mm) were prepared by adding 
3096 H202 and 10% НСІ to remove organic matter and carbonate, 0.5 mol/L Naz2C204 and 0.5 
mol/L (1/63аРОз) were respectively used as the dispersant for neutral and alkaline soils, and the 
sample solution was ultrasonicated for 10 min before transferred to the LS230 analyzer. The 
CaCO3 equivalent was analyzed volumetrically, i.e., by measuring the volume of CO» that was 
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generated by the action of 10% НСІ on soil carbonate. SOC content was determined by the 
K2Cr207-H2S0;4 wet oxidation method of Walkey-Black (Nelson and Sommers, 1982). 

Since Sr and Ca have similar geochemical behaviors, in order to identify the sources of SIC, we 
measured the *"Sr/*6Sr ratio in soil carbonate because it can be used as a proxy for calculating the 
contribution of pedogenic inorganic C (PIC) to total SIC. The method of Chen et al. (1997) was 
adopted to dissolve carbonate in the selected loess and siltigation samples, due to that 0.5 mol/L 
acetic acid solution dissolves carbonate only and does not attack other minerals. In detail, about 
500 mg of 74 um (200 mesh) sample was weighted and dipped in a Teflon container with 20 mL 
of 0.5 mol/L HOAc solution at room temperature for 8 h. Then, the HOAc leachates were 
evaporated in Teflon tubes after addition of 3 mol/L НМОз. The Sr isotopic ratios were 
determined with a Finnigan MAT262 mass spectrometer (Bremen, Germany), installed at the 
Laboratory for Radiogenic Isotope Geochemistry, University of Science and Technology of China. 
All *7Sr/*6Sr ratio values were normalized to the National Institute of Science and Technology 
(NIST) Strontium Carbonate Isotopic Standard (NBS987) value of 0.710288 (30.000013) (26 
standard deviation, n=8). 


2.4 Data collection and statistical analysis 


Climate data were obtained from the dataset produced by the Chinese Academy of Agricultural 
Sciences (http://cdc.nmic.cn/home.do), with a spatial resolution of 1000 m. Data of MAP and 
AMT for each sampling site were extracted on the basis of the geographical locations. 

We measured the *7Sr/96Sr ratio іп soil carbonate as a proxy for calculating the percentage of 
PIC in total SIC (Chen et al., 1997). 


(Us Sr). -0.7080. 
0.004 
PIG =0.12x CaCO,, x PIC; 6, (2) 


where PIC;% is the percentage of pedogenic carbonate in total carbonate of the i horizon; 
(S Sr/*6Sr); denotes the Sr isotopic composition of the i horizon in soil carbonate fraction of the 
0.5 mol/L HOAc dissolved component; PIC; and CaCO;; are the pedogenic inorganic C content 
and CaCO; equivalent content of the 7" horizon (g/kg), respectively; and 0.12 is the conversion 
factor from CaCO; content to inorganic C content. 

This study analyzed the spatial distribution of soil C storage of the top meter. Also, standard 
depths with ranges of 0—30 and 0—50 cm were used to facilitate comparisons with previous studies. 
For each sampling site, the values of SOC and SIC density (kg C/m?) on a volume basis were 
calculated using the following equations: 


PIC% = 100%, (1) 


k 
SOCD , = У p xSOC, x D, x (1- S;), (3) 


i-l 


k 
SICD, = У 0.12х p; x CaCO, x D; x (1- S;), (4) 
i=l 
where SOCD,; апа SICD, are the total amount of soil organic C density and soil inorganic С 
density between soil surface and soil depth h (cm) per unit area (kg C/m?), respectively; k is the 
number of soil genetic horizons within the depth л; p; is the bulk density of the 7" horizon (g/cm?); 
D; is the thickness of the 7" horizon (m); S; is the volume fraction of fragments >2 mm of the i" 
horizon (96); SOC; and СаСОз; are the SOC content and СаСОз equivalent content of the ;'^ 
horizon (g/kg), respectively; and 0.12 is the conversion factor from СаСОз content to inorganic C 
content. 

Multiple comparisons and analyses of variance (one-way ANOVA with Duncan's multiple range 
test) were used to determine the significance of differences between soil C storage among different 
landscapes and among different intervals of soil depth in each landscape. Linear regression and 
correlation analyses were conducted to evaluate the relationship between C storage and climate 
factors (1.e., MAP, AMT). The analyses were performed with SPSS Statistics 20.0 for Windows. 
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3 Results 


3.1 Basic soil properties in different landscapes 


Soils developed in different landscapes in the HRB were classified (Table 1) according to the 
World Reference Base for Soil Resources (IUSS Working Group, 2014), covering most soil types 
in the arid and semi-arid regions of Northwest China. Considerable changes in soil properties 
were observed in different landscapes (Table 2). The generalized genetic horizon sequence of 
soils across the representative landscapes is presented in Figure 2. 


Table2 Basic soil properties of the 10 representative landscapes in the HRB 


Landscape — Horizon Thickness Gravel Clay Silt Sand pH BD SOM CaCO; 
(cm) (%) (%) (%) (%) (g/cm”) (g/kg) (g/kg) 
AM о 12-20 8-10 9-12 55-63 25-36 6.8-7.3 0.88 80.5—82.9 1-2 
A 10-12 15-40 7-13 46-64 23-47 6.6-7.0 0.97-1.11 36.3-67.0 1-2 
B 20-27 10-30 7-8 34-55 38-59 6.4-7.6 1.22-1.33 30.0–39.3 1-2 
2C 57-16 80-85 6-11 26-55 34-68 6.3-7.8 1.15-1.42 7.1-51.1 1-24 
SSM A 14–20 3–8 10-15 48-70 16—42 6.0-8.2 0.72-1.03 . 37.8-127.3 1-16 
B 21–36 5-8 11-15 51-71 14-38 6.2-8.2  0.84-1.17 35.9–85.3 1–30 
2C 26—65 55-85 6-14 37-54 34-57 6.9-8.6 1.32-1.44 4.8–10.9 1-17 
MG A 14-17 0 13-14 53-72 14-34 7.5-8.5  0.87-1.16 56.7–86.3 7-10 
AB 17-28 0 13-15 59-71 16-277 7.4-8.2  1.05-1.17 45.9-59.7 9-61 
B 38-40 0-15 10-16 56-71 14-30 8.2-9.1 1.08-1.22 9.5-47.1 85-192 
C 30-40 0-30 10-14 56-64 27-32 8.8-9.6 1.25-1.32 5.3-7.3 139-153 
MF о 6-15 0 14-19 66-73 13-15 6.9-7.6 0.40-0.50  176.5-203.6 8-11 
A 11-22 0 14-17 65-72 12-20 7.2-7.7 0.40-0.75  115.5–227.0 7-13 
B 32-60 0-2 11-19 65-72 13-23 7.3-7.9 0.57-0.72  132.9-199.0 6-11 
C 10-21 0—50 14-15 69-70 14-16  7.4—7.9  0.61-0.90 | 89.5-178.4 7-14 
2C 13-19 20-50 12-13 63-69 18-25 7.8-8.0 0.88-1.23 28.2-84.8 8-41 
TS A 15-16 2-5 10-12 61-76 128 7.8-8.4 0.91-1.04 42.1—57.5 7-100 
B 66-70 0-10 12-14 64-76 11-23 8.1-8.7 1.07-1.45 5.7–32.5 28-202 
C 25-38 0-10 11-13 67-74 15-20 8.8-8.9 1.23-1.33 7.7-8.1 124-193 
DS A 10-20 0-3 8-12 56-72 16-36 8.3-8.6 1.23-1.26 10.0-11.7 125-139 
B 81-95 0-12 7-13 64-76 15—29 8.0-8.6 1.20-1.32 5.5-9.4 124—160 
С 15-22 0-1 8-10 62-73 17-31 8.2-8.8  1.32-1.36 4.7–7.4 130–143 
HCO A 16-25 1-2 8-15 44-64 21-48 7.9-8.7 1.34-1.45 12.8–23.7 79-105 
B 55-73 1-5 8-16 32-67 21-60  8.0-8.6  1.25-1.55 8.9-17.0 64—147 
С 25—30 1-5 13-14 50-68 18-36 8.1-8.4  1.33-1.45 5.0-11.4 100—124 
RDD A 5-8 10—40 2—14 9-29 57-89 7.6-8.6 1.34-1.46 2.4—5.6 40-76 
C 102-115 2-30 1-6 5-20 74-94 7.5-9.2 1.37-1.57 1.4-2.9 34-56 
AGD A 6-10 1-10 10-12 23-35 53-67 7.9-8.0 1.51-1.67 2.4-3.7 41-100 
B 35-51 10–50 3-14 13-30 56-83 7.9-9.2 1.51-1.77 1.4-2.5 31-71 
2C 48-75 30-90 1-5 7-14 81-92 8.4-9.8 1.55-1.68 1.3-1.6 8-40 
SD С 120 0-2 1-3 4-13 84-95 8.5-9.3 1.26-1.58 0.8–1.2 16–78 


Note: BD, bulk density; SOM, soil organic matter; О, organic horizon; A, surface horizon; AB, transitional horizon; B, subsoil; C, 
substratum; 2C, a lithologic discontinuity with the overlying solum. 


Soils were neutral to strongly alkaline (6.0<pH<9.2) and the fine soil was dominated by the silt 
fraction, except for desert soils in RDD, AGD and SD landscapes, which showed the largest 
proportion in the sand fraction (Table 2). Soil BD increased with increasing soil depth and 
showed large variations among different landscapes, with the lowest BD in MF (0.40-1.23 g/cm?) 
and the highest in AGD (1.51—1.77 g/cm?). In contrast, SOM decreased with increasing soil depth, 
and the highest and lowest SOM contents occurred on MF (28.2—227.0 g/kg) and SD (0.8–1.2 
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g/kg), respectively. Consequently, BD generally showed a significant negative relationship with 
SOM content (r2 —0.91, P«0.001). In addition, the highest and lowest СаСОз contents were 
observed in DS (124—160 g/kg) and AM (1—24 g/kg), respectively. Our results demonstrated 
higher rates of SOM accumulation (i.e., Oi (mattic epipedon and litter layer in AM and МЕ, 
respectively) and Ah (mollic epipedon); Fig. 2) and СаСОз leaching in the wetter and cooler sites 
as compared with those in the drier and hotter sites where the profile includes a calcic horizon 
(i.e., Bk (calcic horizon in the subsoil); Fig. 2) in the HRB. 
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Fig. 2 Generalized genetic horizon sequence of soils in the top meter across the 10 representative landscapes. 
AM, alpine meadow; SSM, subalpine shrub and meadow; MG, mountain grassland; MF, mountain forest; TS, 
typical steppe; DS, desert steppe; HCO, Hexi Corridor oases cropland; RDD, Ruoshui River delta desert; AGD, 
Alxa Gobi desert; SD, sandy desert. Horizon designations follows the methods of Schoeneberger et al. (2012). O, 
organic horizon; А, surface horizon; AB, transitional horizon; B, subsoil; C, substratum; 2C, a lithologic 
discontinuity with the overlying solum. More specifically, Oi, mattic epipedon and litter layer in AM and MF, 
respectively; Ah, mollic epipedon; Ap, anthropic epipedon; Ak, calcic horizon in the topsoil; Bk, calcic horizon in 
the subsoil; Bh, horizon of accumulation of illuvial organic matter; Bu, siltigic horizon. 


3.2 Vertical distribution patterns of soil organic and inorganic carbon in different 
landscapes 


Тће vertical distribution patterns of SOC and SIC density at 20 cm intervals in the top meter are 
presented in Figure 3. Generally, SOC density in all landscapes decreased with increasing soil 
depth (Fig. 3a). There were significant differences of SOC density in the top 20 cm layer and in 
the deeper layers (P«0.05) in different landscapes with exceptions of RDD and SD, where SOC 
density was extremely low (Fig. 3a). The proportion of SOC density in the top 20 cm layer was 
the highest in AM (57% of the SOC density in the top meter), followed by SSM (46%), TS (38%), 
AGD (38%), and MG (35%), suggesting accumulations of SOC in the upper layer in these 
landscapes. In contrast, a much lower SOC density in the top 20 cm layer was found in HCO 
(27%), MF (25%), DS (25%), RDD (23%), and SD (22%), indicating either the importance of 
SOC in the deeper soil layers (1.е., HCO, MF, and DS) or low inputs of SOC in the upper layer 
(i.e., RDD and SD). 

SIC showed different distribution patterns across different landscapes (Fig. 3b). The patterns 
can be categorized into three groups: (1) the relatively uniform distribution of SIC density in soil 
profiles found in AM, SSM, MF, RDD and SD, which exhibited no significant differences of SIC 
density among different layers (P>0.05); (2) the enrichment of SIC in the subsoil of MG, TS, DS, 
and HCO, which had significant higher SIC density in the deeper layers than the upper layer 
(P<0.05); and (3) the accumulation of SIC in the upper layer in AGD, which showed a significant 
higher SIC density in the top 20 cm layer than the other layers (P<0.05). 
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Fig.3 Vertical distribution of (a) soil organic carbon (SOC) density and (b) soil inorganic carbon (SIC) density 
in the 10 representative landscapes of the HRB. Bars mean standard errors (п=3). Different lowercase letters 
indicate significant differences at P«0.05 level among five depth intervals in each landscape. 


3.3 Soil organic and inorganic carbon storage in different landscapes 


Table 3 shows the mean values and standard deviations of SOC and SIC storage at the soil depths 
of 0—30, 0-50, and 0–100 cm. There was a considerable variability in SOC storage and SIC 
storage at each soil depth across different landscapes (Table 3). SOC storage in the top meter 
ranged from 0.82 kg C/m? in SD to 50.48 kg C/m? in MF. The desert soils contained a low SOC 
storage (less than 0.5 kg C/m?) in the top 30 cm soil layer, while MF in forest stored up to 36 
times that amount. The SIC storage at each soil depth (i.e., 0—30, 0—50, and 0—100 cm) was much 
higher in DS than in other landscapes (P«0.05; Table 3). The SIC storage varied from 0.05 to 5.92 
kg C/m? in the 0-30 cm layer, from 0.08 to 10.14 kg C/m? in the 0-50 cm layer, and from 0.19 to 
21.91 kg С/т2 in the 0-100 cm layer. The ratio of SIC storage to STC (soil total carbon) storage 
generally increased with an increase in soil depth (Table 3). The STC storage in the top meter in 
MF (51.20 kg C/m?) was much higher than those («9 kg C/m?) in the barren desert soils (e.g., 
RDD, AGD, and SD). The contribution of SIC to STC varied significantly across different 
landscapes (Fig. 4). The percentage of SIC storage in STC storage in AM, SSM, MG, and MF 
were much lower than those in TS, HCO, DS, RDD, AGD and SD (Fig. 4). 


Table3 Soil organic carbon (SOC) storage and soil inorganic carbon (SIC) storage and the ratio of SIC storage 
to soil total carbon (STC) storage at different soil depths in the 10 representative landscapes of the HRB 


SOC storage (kg C/m?) SIC storage (kg C/m?) SIC/STC ratio 


Landscape 
0-30 cm 0–50 ст 0–100 ст 0–30 ст 0–50 ст 0–100 ст 0-30 ст  0-50cm 0-100 ст 


AM 7.94x1.98* 10.34+2.91°% 11.76+3.141 0.05+0.02° 0.08+0.04° 0.19+0.19° 0.01+0.00 0.01+0.00 0.02+0.01 
SSM 10.89+4.75% 14.42+3.62°° 16.66+0.88° 0.17+0.24° 0.400.62° 0.81+1.30° 0.03+0.04 0.03+0.05 0.04+0.07 
MG 11.33+1.33° 16.52+3.46° 22.37+3.42° 0.92+0.57% 2.86+2.191 12.84+2.32° 0.08+0.04 0.15+0.13 0.36+0.05 
МЕ 18.69+3.80* 29.46+4.30* 50.48+5.35* 0.19+0.07° 0.30+0.09° 0.72=0.22° 0.01+0.00 0.01+0.00 0.01+0.00 

TS 7.68+0.79°% 10.79+0.29° 14.73+1.14°¢ 2.12+1.26°% 5.61+2.97% 15.41+2.11" 0.21+0.11 0.33+0.12 0.51+0.05 
DS 2.06+0.02% 3.24+0.10° 5.83+0.59° 5.92+0.62* 10.14+1.10* 21.9141.61* 0.74+0.02 0.76+0.02 0.79+0.03 
HCO 4.30+0.85% 6.58+1.03% 11.10+1.93* 4.56+0.68° 7.72+0.89° 16.12+0.91° 0.52+0.09 0.54+0.07 0.59+0.06 
RDD 0.50+0.21° 0.78+0.31' 1.54+0.45' 2.0140.25* 3.2630.34* 7.13+0.591 0.80+0.08 0.81+0.07 0.82+0.05 
AGD 0.5140.16° 0.70-0.16° 0.95+0.06' 2.81+1.14° 3.74+1.11% 4.65+0.33" 0.84+0.03 0.84+0.02 0.83+0.02 
SD 0.25+0.01° 0.42+0.03/ — 0.8230.06' 1.88+0.80° 3.2141.13° 6.74+2.27* 0.87+0.05 0.88+0.04 0.88+0.04 


Note: Mean+SE (n=3). Values followed by different lowercase letters within a column are significantly different at P<0.05 level among 
different landscapes for the same soil depth. STC storage was the sum of SOC storage and SIC storage. 
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Fig. 4 Percentage of soil C storage in soil total carbon (STC) storage at the 0—100 cm soil depth in different 
landscapes. STC storage was the sum of SOC storage and SIC storage. 


3.4 Sr isotopic ratio in soil carbonate and its implications for PIC to SIC ratio 


The *’Sr/*°Sr ratio in carbonate through the soil profiles showed a wider ranges from 0.710775 to 
0.711762 and from 0.710965 to 0.711472 in MG and TS, respectively (Fig. 5). In comparison, a 
small variations of *’Sr/*°Sr ratio in carbonate between genetic horizons were found in DS 
(0.711063+0.000059) and HCO (0.711982+0.000038) along their soil profiles (Fig. 5). Following 
the method of Chen et al. (1997), the calculated PIC stocks in the top meter of soils by Equations 
2 and 4 were 7.23, 10.74, 15.50 and 14.59 kg C/m? in MG, TS, DS, and HCO, respectively. 
Correspondingly, the contributions of PIC storage to SIC storage were 7196, 8096, 7696, and 
nearly 100% in these landscapes, respectively, demonstrating the significant role of pedogenic 
carbonate in the formation of SIC in the arid and semi-arid regions irrespective of the limited 
precipitation. 
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Fig.5 °’Sr/®°Sr ratio in soil carbonate along the selected soil profiles in MG, TS, DS and HCO 


3.5 Relationships between soil carbon storage and climatic factors 


SOC storage showed a positive correlation with MAP at different soil depths (1.е., 0—30, 0—50, 
and 0-100 cm), but the coefficient of correlation decreased with increasing soil depth, suggesting 
a diminishing effect of MAP on SOC storage in the deeper layers as compared with the upper 
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layer (Table 4). There were negative correlations between SOC storage (0—30 and 0—50 cm) and 
AMT, and the coefficient of correlation also decreased with increasing soil depth (Table 4). In 
contrast, SIC storage showed an increasing trend with increasing AMT and decreasing MAP. 
However, there were no significant correlations between SIC storage and climatic factors, except 
a weak correlation between SIC storage at the 0—30 cm depth and AMT (Table 4). 


Table 4 Correlations of SOC storage and SIC storage with mean annual precipitation (MAP) and annual mean 
temperature (AMT) at different soil depths 


SOC storage SIC storage 
0-30 cm 0–50 cm 0–100 cm 0–30 cm 0—50 cm 0–100 cm 
MAP 0.82" 0.807 0.68* —0.61 –0.47 –0.22 
АМТ –0.75" –0.68" –0.54 0.64* 0.53 0.35 


Note: * and “* mean significance at P<0.05 and P«0.01 levels, respectively. 


4 Discussion 


4.1 Pedogenic interpretations for contrasting vertical distribution of soil organic and 
inorganic carbon in different landscapes 


The SOC density exhibited a decreasing trend with increasing soil depth in all landscapes in the 
HRB (Fig. 3a). The proportion of SOC density in the 0—20 cm layer relative to that in the top 
meter was the highest in AM (57%) of the HRB, which is similar to the shallow SOC distribution 
in Tibetan alpine meadows (Table 5). This can be attributed to the dense roots in surface soils in 
the alpine meadow (Yang et al., 2010a). It is well known that vegetation, through patterns of 
above- and below-ground allocation, has an important effect on the vertical distribution of SOC 
(Jobbágy and Jackson, 2000). The percentage of SOC density in the top 20 cm layer in SSM (46% 
of the total SOC density in the top meter) of the HRB was much higher than that reported in 
global shrublands (33%; Table 5; Jobbágy and Jackson, 2000). This could be due to the root 
allocation in the shallow layers in subalpine ecosystems and shallow solum with limited addition 
of aeolian deposits (Fig. 2). In contrast, the percentage of SOC in the top 20 cm layer in MF (2596) 
of the HRB was lower than that in global forests (5096; Table 5) (Jobbágy and Jackson, 2000). 
The relatively high proportion of SOC in the deeper soil layers (20-100 cm) in MF could be 
caused by cryoturbation, which leads to the movement of SOM to the subsoils (Baumann et al., 
2009; Liu et al., 2012). The percentage of SOC in the top 20 cm layer in MG (35%) of the HRB 
was lower than that in global grasslands (42%; Table 5). In the semi-arid steppe, the SOC density 
gradually decreased with increasing soil depth (Fig. 3a), which may be related to the deep root 
distribution in this zone. Compared to the average proportion of SOC in the subsoils of the 
croplands in China (63%; Wang et al., 2004), the HCO of the HRB sequestered more SOC in the 
subsoils (7396). This may be related to the deep incorporation of farmyard manure and crop 
residues during the progressive build-up of oases soils. Additionally, the leaching of organic 
matter and colloids along with irrigation has been found to move SOC from surface soils to the 
lower horizons in the drylands (Lal, 2004c). The lack of significant differences in SOC density at 
different soil depths in RDD and SD (Fig. 3a) reflects the homogeneous distribution of SOC in 
soil profiles in the extremely arid environments. However, the relatively high proportion of SOC 
density in the 0—20 cm layer relative to that in the top meter occurred in AGD (38%) highlights 
the physical protection of fine particles in the surface soils of Gobi, which resulted in the 
resistance of organic matter to decompositions in these soils, in comparison with the situation in 
the relatively coarse subsoils (Fig. 2; Table 2). 

The vertical distribution of SIC is in sharp contrast with the decreasing trend of SOC with 
increasing soil depth (Fig. 3). Specifically, AM, SSM, and MF are found mainly in high-altitude 
ecosystems, where intensive carbonate leaching is supposed to occur due to the relatively 
abundant precipitation and low temperature. Furthermore, the higher SOM content in these 
landscapes increases soil porosity through lowering the soil BD (Yang et al., 2014). This could in 
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Table 5 Comparisons of vertical distribution and storage of soil organic and inorganic C in the HRB with 
previous estimates in other studies 
Vertical distribution Soil C storage in the top meter 


Type of soil C (kg C/m’) Reference 
SOC* SIC? SOC SIC 
Grassland 
Global 42% - 11.7–13.2 - Jobbágy and Jackson (2000) 


Ni (2002); Wang et al. (2004); Yu et al. (2007); 


China 39% B 92182 - Mi et al. (2008) 
AM (HRB) 57% (С 11.8 0.2 This study 
Tibetan alpine meadow 5596 - 11.9-12.7, 18.2 11.7, 13.5 Ni (2002); Yang et al. (2010a, b, c) 
Alpine meadow 7 
(Qinghai-Tibet Plateau) - - 25.6 (0-50 cm), 9.5 - Wu et al. (2003a); Chen et al. (2016) 
MG (HRB) 3596 B 22.4 12.8 This study 
Mountain grassland 9.5—20.7 (0-50 cm) И 
(Qilian Mountains) | Y 9.7–19.5 (0–60 cm) 7 ОШ Ziu et el (B007) 
Steppe (HRB) 25%-38% A 5.8–14.7 15.4—21.9 This study 
Tibetan alpine steppe 4196 - 7.4—14.9 14.8 Yang et al. (2010a, b) 
Shrubland 
Global 33% - 8.9 - Jobbagy and Jackson (2000) 
China 46% A 115 - Wang et al. (2004); Yu et al. (2007); Mi et al. 
(2008) 
SSM (HRB) 46% C 16.7 0.8 This study 
Mountain shrub 
(Qilian Novaes) - - 25.1–30.5 (0-50 cm) - Chen et al. (2016) 
Forestland 
Global 50% - 9.3–18.6 - Jobbágy and Jackson (2000) 
China 4296 A 143 - Yu et al. (2007); Yang et al. (2007); Mi et al. 
(2008) 
MF (HRB) 2596 C 50.5 0.7 This study 
Tibetan montane forest 4196 - - - Chang et al. (2015) 
Alpine forest 18.8—31.1 (0-50 cm) ^ 
(Qilian Mountains) ; NGL vue ио 
Farmland 
Global 41% - 11.2 - Jobbagy and Jackson (2000) 
China 34% A 92 Yang et al. (2007); Yu et al. (2007); Mi et al. 
(2008) 
HCO (HRB) 27% A 11.1 16.1 This study 
pana а : : : 16.6-17.1 Wu et al. (2009) 
(Xinjiang) 
Desert 
Global 33% - 6.2 - Jobbágy and Jackson (2000) 
: Yang et al. (2007); Yu et al. (2007); Mi et al. 
China 3296 B 218) - (2008) 
Deserts (HRB) 22960-38706 А, р 0.8-1.5 4.7-1.1 This study 
Gobi (Northwest China) - - 0.7 (0-40 cm) - Zhang and Shao (2014) 
Desert (Hexi Corridor) 3% - - - Wang et al. (2014) 
Desert (Inner Mongolia) - - 9,7 0.1 Wang et al. (2013) 


Note: *, vertical distribution pattern of SOC was expressed by the proportion of SOC density in the 0—20 cm layer relative to that in the 
top meter; >, A represents that SIC content increases with increasing depth, B represents that SIC content peaks at intermediate depths, C 
represents that SIC content shows mixed patterns with increasing depth, and D represents that SIC content decreases with increasing 
depth. -, no data. 


part favor carbonate leaching in forest soils. As a consequence, low SIC density values («1 kg 
C/m?) throughout the soil profiles were found in these landscapes (Fig. 3b). The enrichment of 
SIC in the deeper soils of grasslands (e.g., MG and TS; Fig. 3b) was attributed to the dissolution 
and leaching of carbonates from their topsoil and the subsequent precipitation in the subsoil. This 
is in accordance with the SIC distribution patterns of grasslands in China reported by Mi et al. 
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(2008). However, relatively less carbonate leaching has occurred in the soil profile of DS due to 
the low precipitation (MAP, 144 mm). The homogeneous distribution of SIC in the 0—80 cm soil 
layers in HCO may be associated with the continuous supply of Ca?* and/or Mg?* by irrigation 
and limited carbonate leaching in developing agricultural oases (Gong et al., 1999). Since barely 
any percolating water is available for the migration of SIC in the hyper-arid desert ecosystems, 
parent material might play a major role in shaping the vertical distribution of SIC density for 
RDD, SD, and AGD. The soils of RDD and SD were respectively derived from uniform fine 
alluvium and aeolian sands (Fig. 2), leading to the homogeneous distribution of SIC in both 
landscapes. However, the soils of AGD, which were characterized by fine aeolian deposit rich in 
carbonates overlying the coarse alluvium (Table 1; Fig. 2), exhibited a high percentage of SIC in 
the top 20 cm layer relative to that in the top meter (43%; Fig. 3b). 


4.2 Controlling factors for soil organic and inorganic carbon storage in different 
landscapes 


There was a considerable spatial variability in SOC storage across different landscapes in the 
HRB. Higher SOC storage was observed in high-altitude ecosystems due to the relatively high 
precipitation and biomass production. The highest SOC storage was observed in Picea crassifolia 
forest (MF), which stored more than 50 kg C/m? in the top meter (Table 3). The SOC storage of 
MF at the 0-60 cm depth was comparable to values reported under the same plant species and 
similar climate conditions (Zhu et al., 2017). In this study, we found a relatively low SOC storage 
(11.76 kg C/m?) in AM, which was much lower than that of in Xishui Forest Reserve (Table 5) in 
the Qilian Mountains but was higher than the average of alpine meadow (8.70 kg C/m?) in the 
northeastern margin of the Tibetan Plateau (Chen et al., 2016). These differences could be an 
effect of the considerable spatial heterogeneity across the Qilian Mountains. Although the 
vegetation type was uniform, climate conditions and soil types were very different, thus playing 
an important role in determining the SOC stocks (Liu et al., 2012). The desert soils (RDD, AGD, 
and SD) in the HRB showed the lower SOC storage across the climate transect due to the 
extremely arid climate and low biotic activity. As shown in Table 5, the SOC storage in desert 
soils in our study was similar to that reported in Gobi desert soils (Table 5) by Zhang and Shao 
(2014), but was much lower than the average SOC storage (Table 5) for global desert soils 
(Jobbágy and Jackson, 2000). In the HRB, higher SOC storage was found in HCO as compared 
with that in desert soils, which was attributed to irrigation and fertilization, both enhancing 
biomass production. Overall, the SOC storage in different landscapes generally increased with 
increasing precipitation and decreased with increasing temperature, as evidenced by the 
correlation analyses (Table 4). In drought-prone ecosystems, even a small amount increase of 
precipitation can significantly stimulate plant production and thus contribute to the accumulation 
of SOC (Wu et al., 2003b). Wang et al. (2015) also demonstrated that precipitation is the limiting 
factor of SOC accumulation in the arid and semi-arid regions of Northwest China. The high 
spatial variability of SOC storage in different landscapes in the HRB had important implications 
for ecosystem function and security due to their different sensitivities to climate change. Our 
results also highlighted that an accurate estimation of SOC storage in the arid and semi-arid 
regions of Northwest China should consider the large variations of SOC storage in the diverse 
ecosystems under different landscapes. 

Carbonates have generally been regarded as the most common form of soil C in the arid and 
semi-arid regions (Lal, 2004c; Mi et al., 2008; Monger, 2014), yet little information is available 
on SIC distribution in the HRB. The SIC storage also exhibited large variations across different 
landscapes in this basin. The highest SIC storage occurred in DS (Table 3), due to the arid 
environment together with aeolian deposit rich in СаСОз (Wu et al., 2009). In contrast, the high 
mountain areas (e.g., AM, SSM, and MF) showed a very low SIC storage, which can be attributed 
to the carbonate leaching facilitated by the relative high precipitation. In addition, relatively high 
SIC storage (16.12 kg C/m?) in the top meter was found in HCO due to build-up of siltigation rich 
in carbonate by means of long history of irrigation (Gong et al., 1999). This storage was similar to 
that of irrigated desert soils (16.6 kg C/m?) and irrigated silting soils (17.1 kg C/m?) in China (Wu 
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et al., 2009). Overall, there were no significant relationships between SIC storage and climatic 
factors (Table 4). Previous studies have demonstrated that SIC storage is mainly determined by 
soil parent material and climate (Mi et al., 2008; Wu et al., 2009). The large variations of SIC in 
our study also reflect the differences in parent material and the varying degree of leaching and/or 
calcification of soils. The SIC storage in the arid regions with precipitation less than 200 mm was 
highest in loess (DS), followed by siltigation (HCO), aeolian deposition (surface soil of AGD), 
alluvium (RDD), and aeolian sand (SD; Fig. 6). This indicates that soil parent material is the 
primary control of SIC storage in the hyper-arid environments due to limited leaching of 
carbonate (Wu et al, 2009). Moreover, SIC storage at each depth interval was positively 
correlated with AMT and negatively correlated with MAP (Fig. 7) in the Qilian Mountains where 
soils were derived from aeolian deposits, suggesting climatic factors are important drivers of SIC 
variations in these landscapes (AM, SSM, MG, МЕ TS, and DS). Furthermore, the estimated PIC 
stocks of MG, TS, DS and HCO accounted for more than 70% of total SIC. The highest 
proportion of PIC relative to SIC occurred in HCO due to the abundant irrigation water that 
favored dissolution, translocation and precipitation of CaCOs. The PIC accumulation rate in the 
top meter is 0.6, 0.9, and 1.3 g C/(m?.a) in MG, TS, and DS, respectively, if it is assumed that all 
the carbonates are derived from the Holocene loess (Küster et al., 2006). The low РІС 
accumulation rate in our study agreed well with the results concluded by Scharpenseel et al. (2000) 
and Landi et al. (2003), who reported that the PIC accumulation rate was less than 3 g С/(т2-а) in 
the arid and semi-arid regions of Canada, USA and New Zealand. 

In general, the contrasting landscapes in the HRB lead to greatly different distribution patterns 
of SOC and SIC storage, and there is a general trade-off between SOC and SIC stocks. This 
compensatory effect of environmental conditions, especially climate, on SOC and SIC formation 
can be largely attributed to the fact that the overall decrease in temperature and increase in 
precipitation from arid deserts to alpine mountains facilitate the accumulation of SOC and 
depletion of SIC simultaneously. Further, a better understanding of the counterbalance between 
SOC and SIC would improve the estimation of different C pools and C budgets at watershed 
scales, especially for the arid inland river basins. 
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Fig. 6 SIC storage at three soil depths in regions with mean annual precipitation less than 200 mm. LOE, loess 
(DS); HTM, human-transported material (siltigation; HCO); ALL, alluvium (RDD); EOS, aeolian sand (SD); 
EOD, aeolian deposit (surface soil of AGD). Different lowercase letters indicate statistical significance at Р<0.05 
level among various parent materials for the same soil depth. Bars mean standard deviations. 


5 Conclusions 


We assessed the influences of climate regimes, parent materials and vegetation characteristics on 
SOC and SIC storage and their vertical distribution patterns in 10 representative landscapes in the 
inland Heihe River Basin (HRB), Northwest China. The large observed variations in SOC storage 
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Fig. 7 Relationships between SIC storage at different soil depths with (a) annual mean temperature (AMT) and 
(b) mean annual precipitation (MAP) in the Qilian Mountains where soils were derived from aeolian deposits. * 
means significance at P«0.05 level. 


reflect the very different climatic conditions in different landscapes, which together with the 
corresponding vegetation types, strongly determines the profile distribution. The distribution 
pattern of SIC exhibited more complicated variations than SOC along soil profiles in different 
landscapes in relation to leaching intensity and parent materials at the local sites. Generally, the 
SIC contributed more than 7596 of total C pool when SOC storage was lower than 10 kg C/m? in 
the top meter, suggesting a threshold when SIC cannot be neglected in quantifying C pools in the 
HRB. Furthermore, the *7Sr/?9Sr ratio in soil carbonate proved to be an effective method for 
distinguishing pedogenic and lithologic inorganic C, and PIC accounted for more than 7046 of 
total SIC storage in the selected landscapes, demonstrating the significant role of pedogenic 
carbonate in the arid and semi-arid regions, irrespective of the limited precipitation. The 
considerable variations of SOC and SIC in different landscapes were attributed to pedogenic 
forces resulting from contrasting climatic regimes, parent materials and vegetation types. There is 
a general trade-off pattern between SOC and SIC storage across the contrasting landscapes in the 
HRB due to the changing processes facilitating the accumulation of SOC or SIC in this area. 
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Appendix 


Table S1 Characteristics of sampling sites across the 10 representative landscape types in the Heihe River Basin 


Site Latitude Longitude Altitude (m) Horizon Depth (cm) Gravel (96) 

AMOI 38.62380°N 98.38472°E 4179 Ah 0-10 15 

AB 10-23 15 

Bw 23-43 30 

2С1 43—75 80 

2C2 75-100 80 

AM02 38.228381°N 99.89638°E 3887 Oi 0-12 10 

Ah 12-24 40 

2C 24—100 85 

AM03 38.37852°N 99.32063°E 3630 Oi 0-20 8 

Bw 20-47 10 

2C 47-100 85 

g= SSMOI 38.23413?N 99.89186°E 3682 Ah 0-14 5 

= Ву 14-35 5 

C 2C1 35-57 80 

Г 2C2 57-100 80 

55М02 38.23992?N 99.88897°E 3539 Ah 0-20 5 

Bw 20-42 5 

= 2C1 42-65 75 

2C2 65-100 85 
- SSM03 38.25151°N 99.89244?E 3304 Ah 0-18 
AB 18-38 

Bw 38-74 8 

p: 774—100 55 

MG01 37.89275°N 100.75169°E 3412 Oi 0-16 0 

Ah 16-33 0 

АВ 33-50 0 

x Bk 50-90 15 

а ВС 90—100 30 

b MG02 38.01492°N 100.62564°E 3114 Ah 0-14 0 

É AB 14-42 0 

ВК! 42—63 0 

Bk2 63–80 0 

| Cl 80-113 0 

C2 113-120 0 

MG03 38.02241°N 100.62860°E 3109 Ah 0-15 0 

AB 15-33 0 

Bk 33-60 0 

2Ab 60-90 2) 

2Bkb 90-110 2) 

МЕО1 38.26608°N 99.89329°E 3056 Oi 0-15 0 

Ah 15-37 0 

Bw 37-65 0 

2C 65-78 20 

MF02 38.15239°N 100.25036°E 3004 Oi 0-6 0 

АЋ 6–14 0 

Ah2 14-24 0 

AB 24—40 0 

Bhl 40—53 0 

Bh2 53-75 0 

Bh3 75-100 0 
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Continued 
Site Latitude Longitude Altitude (m) Horizon Depth (cm) Gravel (96) 
MF02 38.15239°N 100.25036°E 3004 Bh4 100-110 0 
Cjj 110-120 0 
MF03 38.11737?N 100.29309°E 2989 Oi 0-11 0 
Ah 11-22 0 
Bwl 22-36 0 
Bw2 36-60 2) 
C 60-81 2 
2C 81—100 50 
TSOL 38.27044°N 99.88842°E 2973 Ah 0-16 2 
ВК! 16–38 0 
Bk2 38-58 0 
Bk3 58—82 0 
C 82-120 0 
TS02 38.22806°N 100.04694°E 2833 Ah 0-15 9 
АВ 15-30 5) 
Bkl 30-52 10 
Bk2 52-70 10 
Bk3 770—100 10 
C 100-120 10 
TS03 38.80214?N 99.94350°E 2739 Ah 0-15 2 
AB 15-25 2 
Bwl 25-46 2 
Bw2 46-70 2) 
Вк 70-95 2 
C 95-120 2 
DSO1 39.14629°N 99.35249°E 2177 A 0-20 1 
ВК! 20–36 1 
Bk2 36-60 1 
Bk3 60-105 1 
C 105-120 1 
DS02 39.26926°N 99.06039°E 2060 A 0-10 3 
ВК! 10–48 12 
Bk2 48—70 p 
Bk3 770—105 1 
C 105-120 1 
5503 39.19134?N 99.48514°E 1879 A 0-17 0 
ВК! 17–47 0 
Bk2 47–98 0 
C 98-120 0 
HCOO1 38.86683°N 100.33540°E 1586 Ар1 0–14 1 
Ap2 14-25 1 
Bul 25-50 1 
Bu2 50-80 1 
Cu 80-110 1 
HCO02 38.88606°N 100.37622°E 1547 Ap 0-16 1 
Bul 16-28 1 
Bu2 28-50 1 
Bu3 50-80 5 
Cu 80-105 5 
HCO03 39.49007?N 99.65261°E 1331 Ap 0-17 2 
Bul 17-30 2 
Bu2 30-54 2) 
Виз 54–78 2 
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Continued 
Site Latitude Longitude Altitude (m) Horizon Depth (cm) Gravel (96) 

HCO03 39.49007?N 99.65261°E 1331 Bku 78-90 2 
Cu 90-120 2 
RDDOI 40.16214?N 99.41900°E 1205 A 0-5 40 
Cl 5-26 15 
C2 26–60 5 
(03 60–80 2 
C4 80—100 S 

C5 100-120 
RDD02 41.03163?N 100.42255°E 1058 A 0-8 10 
2С1 8—20 30 
2C2 20–58 30 
2C3 58-100 3 
RDD03 40.77397?N 100.13444°E 1111 A 0-8 10 
C 8—38 15 
T 2C1 38-73 15 
2C2 73-103 p 

AGDOI 41.87041?N 100.58152°E 965 A 0-6 

Bwl 6-30 15 
Bw2 30-57 50 
2С1 57-86 90 
2C2 86-120 90 
AGD02 41.69688°N 100.28950°E 992 Ay 0-10 10 
BC 10-52 30 
2С1 52—82 30 
2C2 82-100 30 
AGDO03 41.23336?N 100.50868°E 1033 A 0-10 3 
Bw 10-45 10 
2С1 45—86 75 
2C2 86-120 75 
50101 39.35686?N 99.14116°E 1600 Cl 0-9 0 
(62 9-85 0 
C3 85-120 0 
= SD02 39.62055°N 99.60005°E 1312 Cl 0-10 2 
$ C2 10-40 2 
C3 40-100 2 
5103 40.66165?N 101.42041°E 1068 Cl 0-15 2 
(62 15-65 2 
(65 65-120 2 


Note: AM, alpine meadow; SSM, subalpine shrub and meadow; MG, mountain grassland; MF, mountain forest; TS, typical steppe; DS, 
desert steppe; HCO, Hexi Corridor oases cropland; RDD, Ruoshui River delta desert; AGD, Alxa Gobi desert; SD, sandy desert. 
Horizon designation follows the methods of Schoeneberger et al. (2012). Horizon: O, organic horizon; A, surface horizon; B, subsoil; C, 
substratum; AB and BC, transitional horizon. Horizon suffix: b, buried genetic horizon; h, organic matter accumulation; i, slightly 
decomposed organic matter; jj, evidence of cryoturbation; К, pedogenic CaCO3 accumulation; p, plow layer; u, presence of 
human-manufactured materials; w, weak color or structure within B (used only with B); y, accumulation of gypsum. Numerical prefix 
(i.e., 2) is used to denote lithologic discontinuity (e.g., 2C, 3C). Numerical suffixes (i.e., 1, 2, 3) are used to denote subdivisions within a 
master horizon (e.g., ВК!, Bk2). 


